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Using inelastic neutron scattering, x-ray, neutron diffraction, and the first-principle calcula- 
tion techniques, we show that the crystal structure of the two-dimensional quantum spin system 
(CuCl)LaNb207 is orthorhombic with Pbam symmetry in which CUCI4O2 octahedra are tilted from 
their high symmetry positions and the Cu'^'*'(s = 1/2) ions form a distorted square lattice. The dom- 
inant magnetic interactions are the fourth nearest neighbor antiferromagnetic interactions with a 
Cu-Cl-Cl-Cu exchange path, which lead to the formation of spin singlets. The two strongest in- 
teractions between the singlets are ferromagnetic, which makes (CuCl)LaNb207 the first system of 
ferromagnetically coupled Shastry-Sutherland quantum spin singlets. 



PACS numbers: Vaiid PACS appear here 

Quantum magnets continue to attract interest due to 
the wide range of exotic ground states and excitations 
that they display, [ij Such magnets usually have mag- 
netic ions with low spin values (typically spin-1/2 or 
1) which are coupled by low-dimensional and/or frus- 
trated magnetic interactions. A quantum magnet that 
has attracted a lot of theoretical interest is the square 
lattice of spin-1/2 ions where the first (Ji) and second 
( J2) neighbor interactions which couple the magnetic ions 
along the sides and diagonals of the square lattice, respec- 
tively, compete with each other. However in many 
cases the predictions have not been verified experimen- 
tally due to the lack of physical realizations of this sys- 
tem. Recently (CuCl)LaNb207 was proposed as an ex- 
ample of frustrated square lattice. Here the magnetic 
Cu^"*" ions which possess spin moment S=l/2 are octahe- 
drally co-ordinated and form CuCl planes well separated 
by LaNb207 layers. The crystal symmetry was originally 
found to be tetragonal with space group P4/mmm and 
one copper ion per unit cell {at = bt = 3.88 A) forming an 
ideal undistorted square lattice of S — 1/2 ions with uni- 
form Ji and J2.[3 Theory suggests that (CuCl)LaNb207 
would have long-range magnetic order which is ferromag- 
netic or antiferromagnetic depending on the ratio of the 
two exchange constants. Q Gapped excitations are not 
expected except for a very limited region of the phase 
diagram when the ratio of interactions along the diago- 



nal to those along the square are w O.sQ while a gap- 
less spin disordered phase is predicted when this ratio 
is w -0.5. Q However the experimental data appears to 
contradict the theory, heat capacity and magnetization 
measurements suggested that the ground state is in fact 
a singlet and that the excitations are gapped. [1| Neutron 
scattering also revealed the absence of long range mag- 
netic order and showed that the excitations are gapped 
and centered at 2.3 meV.Q Such behavior is typical of 
a dimerised magnet where a dominant antiferromagnetic 
exchange interaction couples the spins into dimer pairs 
with no net moment and a finite energy is required to 
break up the dimers and create a triplet (5 = 1) exci- 
tation. However such dimerization is incompatible with 
the high crystal symmetry where Ji and J2 are uniform 
and the spins cannot be paired. 

This paper describes a reinvestigation of the crystal 
structure and shows that the space group is in fact or- 
thorhombic (Pbam) with a unit cell that is doubled along 
both directions within the afe-plane compared to the orig- 
inal tetragonal symmetry, thus allowing for dimerisation. 
The Cu-Cl bonds are considerably distorted and band 
structure calculations reveal highly spatially anisotropic 
exchange interactions and predict dimerisation between 
fourth neighbour ions, in agreement with our inelastic 
neutron scattering data. The resulting magnetic struc- 
ture is best described as antiferromagnetic dimers which 
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FIG. 1: Inelastic neutron scattering data obtained at 70 mK 
from the powder sample, (a) Energy dependence of the triplet 
excitations for B = T (blue) and B = 11.5 T (red), (b) Q- 
dependence for B = T. Lines are described in the text. 



are coupled together by frustrated ferromagnetic interac- 
tions in a Shastry-Sutherland type arrangement. 

A 15 g powder sample of (CuCl)LaNb207 was synthe- 
sized by an ion-exchange reaction using a powder mix- 
ture of RbLaNb207 and CuCl2 with a molar ratio of 1:2 
heated at 335 °C for one week, as described previously [H. 
Concerning single-crystal samples, precursor CsLaNb207 
single crystals (see Ref. jlo|) with typical size of 0.5 x 
0.5 X 0.1 mm'^ were embedded in a molar excess of CuCl2 
powder (99.99%) sealed in an evacuated tube and reacted 
at 340 °C for one week. The mixture was then washed 
with distilled water to remove the byproducts and un- 
reacted CuCl2. The square, plate-like dark green crys- 
tals of (CuCl)LaNb207 were then dried overnight at 120 
°C. The Energy Dispersive Spectroscopy (EDS) study re- 
vealed the absence of the Cs atoms, indicating that the 
ion exchange reaction occurs completely even for the sin- 
gle crystal. 

A CCD and a 4-circle diffractometers were used for 
the single-crystal data collection at 293 and 14 K, re- 
spectively. Diffraction patterns of the powder were mea- 
sured at room and low temperature using synchrotron 
light with A = 0.77749 A at BL02B2 bcamline at SPring- 
8 as well as neutron beam with A = 1.9085 A at DIA 
at ILL. Powder inelastic neutron scattering measure- 
ments were done with A = 3.8 A at the Disk-Chopper- 
Spectrometer (DCS) at NCNR, NIST under an external 
magnetic field. The crystal structure was solved by direct 
methods against single crystal data and further refined 
by taking into account the twinning due to the pseudo- 
symmetry. The structural parameters were optimized by 
a Reitveld refinement against neutron diffraction data, 
more sensitive for O and CI atoms. All the refinements 
were carried out by the means of the JANA2000 refine- 
ment program [ll| . 

Fig. 1 (a) shows the inelastic neutron scattering data 



TABLE I: Structural parameters obtained from Rietveld re- 
finement of neutron diffraction at 2 K. The goodness of the 
fit is indicated by the conventional structural Rsragg, Rwp, 
Rp, and x^- 

Pbam, a = 7.7556(5) A, fe = 7.7507(5) A, c = 11.7142(4) A 
= 3.38, Rflragg = 1.64%, Rp = 3.39 %, R^p = 4.62 % 



Atom(iy) 


x/a 


y/b 


z/c 


f/iso(A^) 


Cu (4ft) 


0.502(6) 


0.7662(10) 


0.5 


0.0043(13) 


CI (4ft) 


0.7367(16) 


0.4286(7) 


0.5 


0.006(2) 


La (4p) 


0.7412(11) 0.5001(10) 





0.0024(12) 


Nb (Si) 


0.500(3) 


0.7529(9) 


0.1908(3) 


0.0040(9) 


Ol (4/) 


0.5 





0.1397(13) 


0.002(3) 


02 (4e) 


0.5 


0.5 


0.1751(11) 


0.008(3) 


03 (Si) 


0.752(4) 


0.7509(11) 0.1545(10) 


0.006(2) 


04 (4^) 


0.500(2) 


0.7170(14) 





0.007(2) 


05 (Si) 


0.501(4) 


0.7731(11) 


0.3410(5) 


0.0089(16) 



as a function of energy transfer, Hu. In the absence of 
an external magnetic field, i? = T, a single peak is ob- 
served centered at Huj = 2.22(1) meV in agreement with 
previous work.^ An external magnetic field of i? = 11.5 
T splits this peak into three peaks at hu = 1.19(1), 
2.45(3), and 3.95(11) meV. The energies are consistent 
with the Zeeman splitting expected for a magnon excita- 
tion with spin S = 1 {g^sSzE) thus confirming that the 
excitations are due to a transition from a singlet ground 
state to a triplet excited state and implying that the 
spins in (CuCl)LaNb2 07 are paired into dimers where 
the dimer exchange constant is « 2.22 meV. In zero field 
the excitation extends in energy from ~ 1.2 meV to ~ 3.0 
meV and is thus much broader than the resolution (0.16 
meV) implying that it is dispersive due to interactions 
between dimers. Indeed the lower edge of the excita- 
tion is consistent with the critical field of Bd = 10.3 T 
from magnetization measurements, for condensation of 
magnons into the ground state. [Sj] 

The nature of the spin dimers in (CuCl)LaNb207 
can be obtained from the wavevector-dependence of the 
triplet excitations, I{Q). As shown in Fig. 1 (b), 
the Q-dependence exhibits a prominent peak centered 
at Qc ~ 0.5 A^^. If we assume the dimers are non- 
interacting, the Q-dependence of their excitation goes 

as, 1 — SlsSl where r is the intradimer distance. The 

' Qr 

best fit of the data to the isolated spin dimer model was 
obtained with r = 8.4(2) A (black line in Fig. 1 (b)) 
in agreement with previous work.fij The model repro- 
duces the observed peak position of Qc ^ 0.5 A^^ well. 
This reveals that the dominant magnetic interactions are 
between the Cu'^^ ions that are separated much further 
than the nearest neighbor distance of 3.88 A, a result 
that is inconsistent with the currently accepted structure 
of (CuCl)LaNb207. 

To understand this unusual Q-dependence, we have 
reinvestigated the crystal structure in great detail. Fig. 2 
(a) shows the CCD image of the diffraction data taken for 
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FIG. 2: Single crystal laboratory X-ray and powder neutron 
diffraction data of (CuCl)LaNb207. (a) the X-ray CCD im- 
age of the afo-plane at 14 K. The single crystal was twinned, 
thus a- and fe-axes are interchangable. The spots with green 
circles are the allowed integer Bragg reflections for the origi- 
nal tetragonal structure, the other spots are the half-integer 
reflections, (b) Neutron diffraction pattern at 2 K. Dots are 
the data, solid line is the fit to the orthorhombic structure, 
and the line on the bottom panel is the difference between 
the data and the model. The indices of the reflections in the 
inset are in the tetragonal notation. 





FIG. 4: (a) Magnetic exchange interactions in the a6-plane 
of (CuCl)LaNb207. The blue, green, and red spheres repre- 
sent the Cu, CI, and O ions, respectively. The O atoms are 
located on the top and bottom of the Cu atoms. The lines 
connecting Cu atoms represent exchange bonds: Jia (blue), 
J2a (solid cyan), J2b (dotted cyan), J4 (solid red), and J4 (dot- 
ted red line). See Table II. The negligible nearest neighbor 
interaction along the 6-axis, Ju, is not plotted for clarity, (b) 
Six ordered antiferromagnetic spin states used to extract the 
spin exchange parameters by GGA+U calculations. Only the 
Cu^^ ions are shown, and the open and filled circles represent 
the up- and down-spins, respectively. The relative energies of 
the ordered spin states are given in the parentheses. 



TABLE II: Exchange paths for up to the 4th nearest neighbor 
interactions in (CuCl)LaNb207 and the corresponding cou- 
pling constants relative to J4. J4 is antiferromagnetic and of 
strength ~ 2.22 meV. Negative ratio Js/Ji means that Js is 
ferromagnetic. Short and long lines between atoms represent 
different bond lengths without proportion. 



Js 


path 


d(A) 


angle (°) 


Js/Ja 


Jla 


Cu-Cl-Cu 


3.88548 


108.9, 75.8 


-0.39 


Jib 


Cu-Cl-Cu 


3.88548 


80.9 


-0.04 


J2a 


Cu-Cl-Cu 


5.46148 


156.7 


-0.38 


J2b 


Cu-Cl-Cu 


5.5053 


170.2 


-0.14 


Ja 


Cu-Cl-Cl-Cu 


8.81262 


164.9 


1 


J'i 


Cu-Cl-Cl-Cu 


8.53250 


150.0 


0.18 



FIG. 3: Crystal struture of (CuCl)LaNb207. (a) Three 
dimensional view, (b) the fee-projection, and (c) the ac- 
projection. Purple, red, blue and green spheres represent La, 
O, Cu and CI atoms, respectively. The dark green polygons 
represent NbOe octahedra. The grey lines show the unit cell. 



the (hkO) plane at 14 K. Note that the crystal is twinned 
and thus h and k are interchangeable. In addition to 
the allowed integer peaks for the PA/mmm tetragonal 
structure, half-integer superlattice spots are also found, 
indicating that the real unit cell is 2at x 2bt x c com- 
pared to the tetragonal unit cell. Similar superreflections 
are also observed in the powder synchrotron x-ray and 
neutron diffraction (see Fig. 2 (b)). Furthermore, the 
CCD image shows that the superlattice peaks (odd hflfl) 
and (O,odd k,0) are absent in the orthorhombic notation 
(red arrows). This extinction condition tells us that the 
crystal structure is orthorhombic with either Pbam or 



Pba2 space group. The best fit to the x-ray and neutron 
diffraction data was obtained with the Pbam structure. 
Table I lists the optimal structural parameters from the 
neutron diffraction data that are more sensitive to the 
CI and O atoms than the x-ray data. The validity of 
the structural analysis is further supported by the bond 
valence sum calculation for copper yielding +2.02. A 
comparison of diffraction patterns at different tempera- 
tures revealed that no structural phase transition occurs 
below 293 K. 

The refined orthorhombic crystal structure differs 
markedly from the tetragonal structure previously re- 
ported. While the NbOe octahedra align with the axes in 
the tetragonal structure they are strongly tilted in the or- 
thorhombic structure particularly around the a-axis in a 
staggered manner. Furthermore, the Nb ions are shifted 
along the c-axis from the center of the octahedron. The 
tilting pattern of the NbOg strongly influences the po- 
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sitions of both the copper and chlorine atoms. In par- 
ticular, the chlorine atoms move significantly along the 
^-direction and slightly along the a-direction from their 
tetragonal position, probably to reduce Coulomb repul- 
sion between chlorine and apical oxygen atoms. The Cu 
ions occupy the Ah sites. Along the a-direction, they 
are mostly transversally displaced, i.e., along the 6-axis, 
yielding the distance between the nearest neighbor Cu- 
Cu ions to be 3.626 A and 4.129 A along the 6-axis and 
3.885 A along the a-axis. The copper ion is coordinated 
octahedrally by two oxygen ligands with a distance of 
1.865 A as well as four chlorine ligands with two shorter 
bonds (2.386 A, 2.389 A) and two longer bonds (3.136 A, 
3.188 A). When the local z- and x-axes for each Cu^'*" ion 
are taken along the Cu-0 and the short Cu-Cl bonds, re- 
spectively, the overall symmetry of the magnetic orbital 



including the ligand p-orbitals has the z 



character. 



As a result, the spin exchange interactions in the CuCl 
layer become highly anisotropic. 

As for the spin exchanges of (CuCl)LaNb207, we con- 
sider Cu-Cl-Cu and Cu-Cl-Cl-Cu superexchanges. In 
general, one expects that the Cu-Cl-Cu exchange is fer- 
romagnetic if the bond angle is close to 90°, and anti- 
ferromagnetic if it deviates from 90°. However, when 
a Cu-Cl-Cu exchange path contains a long Cu-Cl bond 
(represented by Cl-Cu-Cl in Table II), the exchange be- 
comes ferromagnetic even if the Cl-Cu-Cl angle deviates 
considerably from 90° because the magnetic orbital is 
not contained in the long Cu-Cl bond. IJ] The Cu- 
Cl-Cl-Cu spin exchange should become more strongly 
antiferromagnetic with increasing Cu-Cl-Cl angle and 
shortening Cl-Cl contact distance so that the overlap 
between the CI 3p orbitals in the Cl-Cl contact becomes 
large. ^l3| The six exchanges considered in our study are 
listed in Table II and described in Fig. 4 (a). To eval- 
uate these exchanges, we determine the relative energies 
of seven possible ordered spin states of (CuCl)LaNb2 07 
(see Fig. 4 (b)) on the basis of density functional cal- 
culations employing the frozen-core projector augmented 
wave method UA UM encoded in the VASP Il7| with 
the generalized-gradient approximation (GGA) for the 
exchange-correlation functional 18| . The GGA plus on- 
site repulsion U (GGA-I-U) method [13] with effective U 
— 4 eV was used to properly describe the strong elec- 
tron correlation of the Cu 3d states. By mapping the 
relative energies of the seven states determined from the 
GGA-f U calculations onto the corresponding energies de- 
termined from the spin Hamiltonian defined in terms of 
the six exchanges, we obtain their values listed in Table 
II. The fourth nearest neighbor (NN) interaction, J4, of 
the Cu-Cl-Cl-Cu exchange type is the strongest and an- 
tiferromagnetic. The other fourth NN coupling J4 is also 
antiferromagnetic, but is much weaker than J4: J4/J4 = 
0.18. This is because J4 has a larger Cu-Cl-Cl angle and 
a shorter Cl-Cl distance than J^: 164.9° and 3.835 A for 
J4 vs. 156.0° and 4.231 A for j',. AU other Js are fer- 



romagnetic. The strengths of the six exchanges decrease 
in the order, J4 > Jia > J2a > J4 > J2b > Jib- 

The fact that J4 is both the strongest interaction and 
is antiferromagnetc explains the long standing mystery 
of the spin singlet formation in (CuCl)LaNb207; the dis- 
tance between the Cu ions connected by J4 is indeed r = 
8.533 A, consistent with the intradimer distance obtained 
by the isolated spin dimer model. To improve the model 
we considered the exchange paths Jia and J2a in addition 
to J4 and fitted the inelastic neutron scattering data to 
the first moment sum rule for powder data [20| 



In . 



a-^ J,(So-SdJ|/cu^+(g)|' (1- 



sin{Qds 
Qds 



In this equation the integration is over solid angle and en- 
ergy, Js is the exchange constant coupling the s*^ nearest 
neighbor spins, d^ is their separation and (Sq.Sci^) is the 
two spin correlation function for this pair. From this 
expression it is clear that each coupling constant Jg pro- 
duces a modulation in the first moment with a periodicity 
depending on the separation d^ of the spins. Therefore 
the dominant exchange interactions can be deduced but 
the magnitude of the exchange constant cannot be de- 
termined because the fitted quantity is the product of 
the exchange constant and the spin correlation function. 
The best fit was obtained with J4 {So-SdJ = 0.027(2), 
Jia(So.SdiJ = 0.011(2) and J2a{So.Sd^J = 0.005(2). 
Note that the ground state obtained by the first princi- 
ple calculations is the AF3 state shown in Fig. 4 (b), 
and thus a positive value of Js (So.S^^) indicates that 
the spin correlation agrees with the sign of the exchange 
constant, i.e., the spins are parallel for a ferromagnetic 
interaction and antiparallel for an antiferromagnetic in- 
teraction. The fit is much better than that of the simple 
dimer model, and the fitted parameters are consistent 
with our first principle calculation results. This indi- 
cates that in (CuCl)LaNb207 the fourth nearest neigh- 
bor Cu^^ ions form spin singlets, which are arranged 
orthogonally in the a5-plane, while the coupling between 
them is primarily ferromagnetic. Thus, the spin lattice 
of (CuCl)LaNb207 is best described as ferromagnetically 
coupled Shastry-Sutherland quantum spin singlets. 
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